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Abstract 

We perform a time-dependent study of the driven dynamics of overdamped particles which are 
placed in a one-dimensional, piecewise linear random potential. This set-up of spatially quenched 
disorder then exerts a dichotomous varying random force on the particles. We derive the path 
integral representation of the resulting probability density function for the position of the particles 
and transform this quantity of interest into the form of a Fourier integral. In doing so, the evolution 
of the probability density can be investigated analytically for finite times. It is demonstrated that 
the probability density contains both a <5-singular contribution and a regular part. While the 
former part plays a dominant role at short times, the latter rules the behavior at large evolution 
times. The slow approach of the probability density to a limiting Gaussian form as time tends to 
infinity is elucidated in detail. 

PACS numbers: 05.40.-a, 05.10.Gg, 02.50.Ey 
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I. INTRODUCTION 



The Langevin equation for an overdamped particle in a one- dimensional random potential 
presents a useful model for the study of a variety of the noise-induced and disorder-induced 
phenomena. Specifically, it has been employed for describing anomalous diffusion in dis- 



ordered media 



assy dynamics and depinning transition 



transport in ratchets [9, 



mobility in disordered solids 




containing explicitly quenched disorder 



7. 



13 



, diffusive 



14| , carrier 



lit] ], to name only a few main ones. Moreover, this set- 
up enables a very simple description of the dynamics of domain walls in random magnets, 
dislocations in solids, and vortices in type-II superconductors, effectively capturing all the 
essential effects of quenched disorder. 

The effects arising from quenched disorder are often unexpected and even counterintu- 
itive. It therefore constitutes a prominent objective to study model systems for which exact 
statistical characteristics of the solution of the corresponding Langevin equation can be ob- 
tained, not only at asymptotically large times but even at finite times of its evolution. Exact 
results are very useful also for testing approximate methods of analysis of these systems. 
There already exists in the literature a number of exact results derived for systems with 



Sinai disorder 
order 



isorder jj, 2. 17. 

bhhBq 



.7^Q, Q, Gaussian disorder jj], ^, [(j, [g, [20] and non-Gaussian dis- 



In most of these cases the obtained results represent either the 
numerical characteristics of the solution of the Langevin equation (moments of the solution, 
moments of the first-passage time, etc.) or the asymptotic probability distributions describ- 
ing the long-time behavior of the solution. To the best of our knowledge, however, the exact 
time- dependent distributions of analytic form, allowing the study of the effects of quenched 
disorder in full detail, have not been given before. 

In this paper we solve the above-mentioned challenge for a special class of systems contain- 
ing quenched dichotomous disordered and subjected to a constant bias force. More precisely, 
we derive and investigate the analytical expression for the probability density function of 
overdamped particles in a piecewise linear random potential driven by a constant force at 
zero temperature. The paper is organized as follows. In Sec. II, we describe the model and 
introduce the path integral representation of the probability density. In Sec. Ill, we reduce 
the probability density to the form of a Fourier integral. The first two moments of the proba- 
bility density and their asymptotic expressions are calculated in Sec. IV. The time evolution 
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of the probability density is studied both analytically and numerically within Sec. V. Here 
we also compare the analytical results with those obtained from direct numerical simulation 
and derive the asymptotic behavior of the probability density in the limits of small and large 
times. In Sec. VI, we summarize our findings. Some technical details of our calculations are 
deferred to the Appendixes A and B. 

II. BASIC EQUATIONS AND DEFINITIONS 

We consider the overdamped dynamics of a classical particle governed by the dimension- 
less equation of motion 



Here, X t denotes the particle coordinate, /(> 0) is a constant force, and g(x) = —dU(x)/dx 
is the dichotomous random force produced by the piecewise linear random potential U(x) 
(see Fig. 1). We assume that U(x) is characterized by (i) statistically independent random 
intervals of lengths Sj which are distributed with the identical probability density p(s) and 
(ii) two deterministic potential slopes +g and —g, so that g(x) = =fg. In addition, we assume 
that at the initial time t = a particle is located at the origin of the coordinate system and 
that all sample paths of U(x) start with a positive slope dU(x)/dx\ x=+ o = +g, so that 



Under these assumptions the dynamical solution of Eq. (12. ip exists only if / > g; otherwise, 
if / < g, the particle stays localized in its initial well. 

Equation (12. ip is of minimal form in order to account for the effects of quenched disorder 
on the overdamped motion of driven particles. As will be demonstrated below, its main 
advantage is that many of the statistical properties of X t can be derived exactly. Moreover, 
Eq. (12.11) can be used also for studying a number of important physical issues ranging from 
the low-temperature dynamics of charge carriers in randomly layered media to the light 
propagation in interstellar space. 

Our main objective is to express the probability density Pf(x) that X t = x for a fixed t 
via the dimensionless characteristics of the force field. In order to calculate this quantity, 
we start out from its definition 



X t = f + g(X t ). 



(2.1) 



X = 0, g(+0) = -g. 



(2.2) 



P t (x) = (S{x-X t )), 



(2.3) 
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wherein the angular brackets denote an average over the sample paths of g(x) and S(x — X t ) 
is the Dirac 5-function. We first introduce the total probability W n (t) of those sample paths 
which have n(> 1) jumps on the interval (0, X t ). In this case the particle coordinate X t lies 
between the nth and (n + l)st jumps of g(x) and thus can be represented as 



1 1 



X t = J2 S 3 + § n+l (2.4) 
3=1 

where s n+ \ < s n+ \. Because, according to Eq. (12. ip . the particle passes the interval Sj during 
the time Sj/[f + (— lYg], the relation 

§/+(-i)^ + /+(-V +i ^ =t (2,5) 

holds. It implies that 

Si 



The condition 



+1 = [/ + (-l)- 9 ]( t -g 7Tlfi5 -j. (2.6) 



being a consequence of the relation s n+ i > 0, defines in the n-dimensional space of the 
parameters Sj a domain Q n (t) for their allowed values. Therefore, since the probability that 
the {n + l)st jump of g(x) occurs at s > s n+ i is equal to p(s)ds, we readily obtain 

W n (t)= (U ds A^) / P(s)ds. (2.8) 

JQ„(t) \ j—i / JSn+l 

If the function g(x) has no jumps on the interval (0,X t ), i.e., if n = 0, then the total 
probability of these sample paths is given by 

rod 

W (t) = / p(s)ds, (2.9) 

J si 

where sj = (f—g)t. We note that the probabilities W^o(t) and W n (i) are properly normalized: 
W Q (t) + E^=i ^n(t) = 1 (see Appendix A). 

The above results, together with the definition A2.3I) . show that the probability density 
function can be written as 



P t (x) = 5[x - (/ - y)t]W (*) + E p t\x) , (2-10) 
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where 




n 

x8(x-J2 s J -Sn+i)- (2-11) 

3=1 

These formulas provide a path integral representation of Pt(x) in the case of a dichotomous 
varying random force g(x). We emphasize that this representation is quite general and 
remains valid for an arbitrary probability density p(s) of the random intervals Sj. 



III. TIME-DEPENDENT PROBABILITY DENSITY FUNCTION 



From a practical point of view, the result (12.101) is inconvenient because of its com- 
plexity. The main difficulty arises from the necessity to perform the integration over the 
n-dimensional domain fl n (t) and the summation over all n. Here we solve this problem in 
the case when the lengths of the random intervals Sj are distributed with the exponential 
probability density p(s) = Aexp(— As), where A -1 denotes the average length of these inter- 



vals. The key point o: 
unit step function [23] 



our approach consists in the use of the integral representation of the 



1 

2^ 



+oo p (iz+y)K 



iz + y 



dz 



1 if k > 



if k < 



(3.1) 



where y is a positive-valued, real parameter. Applying (13. 1 p to (12. lip and setting k 
t - £?=i Sj/[f + (-l) j g], we obtain 



Pt\x) 







dz . 




/-oo iz + y Jo 





xe 



-(i z +y) Sj /[f+(-l)ig] 



n ds jP ( Sj ) 

3=1 

p(s)ds 



Sn + l 



3=1 



(3.2) 



Though in (13.21) an extra integral appears, this form of (x) is much more advantageous 
as compared to the integration over the domain Q n (t) in ( 12 . 101) because it contains the 
independent integrations over the n variables Sj. Using the integral representation of the 
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*(0 = 4 r e"" 4 * ** (3.3) 



oo 



£-f unction 

am = 

and taking into account that according to ( 12. 5ft 

X) s i + s 2m+i = (/-ff)t + 7— Ss?,, (3.4) 
i=i J + 5 i=i 

the formula (13.21) for p(s) = Aexp(— As) and n = 2m can be rewritten in the form 



00 



Here, 



with Re i/ > and 



Similarly, since 



(27r) 2 J-oo J-oo iz + y 

xe [iz+y ~ x{f - 9 ^ t I m (iy l )I m (iy 2 ). (3.5) 



I(v) = / p(s)e- ( "- A)s rfs = - (3.6) 
Jo V 



iz + y iz +y + 2g(X-irj) 

" 1 = 7^- " 2 = — rr a — ■ (37) 



2m— 1 2^ m 

53 Sj + S 2m +1 = (/ + g)t ~ J— 53 S 2j-h (3-8) 



the formula (13.21) for n = 2m — 1 yields 



(27r) 2 J-oo J-oo iz + y 

xe [tz+y - x{f+9)]t I m (u 3 )I m -\u 4 ) (3.9) 



with 



iz + y - 2g(X - irj) iz + y 

= 7 > ^4 = -7- — • (3.10) 

Using the method of contour integration, it is easy to show that the integrals over z in 
( 13. 5p and (13.91) do not depend on the arbitrary positive parameter y. This means that, in ac- 
cordance with the definition ( 12 . 1 If) . the partial probability densities P^ 2m \x) and pj> 2m ^(x) 
do not depend on y as well. Therefore, for calculating the probability density Pt(x) we can 
use the most appropriate values of this parameter. Specifically, if before the integration over 
z we want to sum over m, it is reasonable to choose y > X(f + g) in order to avoid dealing 
with divergent sums and integrals. In this case the condition Re Vj. > A holds for all k, and 
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thus |I(^fc)| < 1- The last inequality permits one to use a geometric series formula for the 
evaluation of series like 

t /•>,)/•>,) = . y, ( *> , , 

Ei-wi-'w = r^Miwy (3 ' n) 

Then, substituting W (t) = e" A(/ " 9)i and 



oo 



E pr(x) = £ \pr\x) + py^ix)} (3.i2) 



n=l m=l 

into (12.101) and collecting the above derived results, we find the probability density Ptix) in 
the form of the Fourier integral, i.e., 



P(x) = ±- r Mv)e- ir > x d V , (3.13) 

Z7T J-oo 



where 



a -(*-">)/* r°° liz + y + 2g(X - irf) 



2tt 



e {\-irf)gt 



I Kf-9) e -(\-i V )gt 
(z - z 3 )(z - z 4 ) 



(z - z x ){z - z 2 ) 

e (iz+v)t dz (3.14) 



is the characteristic function of X t , 

z\ = gr] + a{rj) + i[y + Xg + b(rj)], 
Z2 = grj- a (rt) + i[y + Xg - b(rj)], 
Z3 = -grj + a(v) +Av~^9 + Hv)], 

z 4 = -gij- a(rj) +i[y - Xg -b(rj)], (3.15) 

and a(rj) = Xg 2 r]/b(r]), 

Kv) = ^{A 2 / 2 - gW + [(A 2 / 2 - gWf + 4AW]*}*. (3.16) 

Finally, calculating the integral in (13.141) by the method of contour integration (see Ap- 
pendix B), we find the characteristic function 4>t{v) i n terms of elementary functions 



( v ) = e -(*-iv)ft ( l 9V sin h[6( 7y )t - ia0n)t] 

\b(r)) - ia{rj) 



+ cosh[b{r])t-ia(r])t]). (3.17) 
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As expected, the characteristic function does not depend on the auxiliary parameter y. 
The functions a(rj) and b(r]) tend to grj and Xg, respectively, as |r/| — > oo. Accordingly, 
4> t (r]) oc eWs* if | 77 1 — > 00, and so the probability density Pt(x) contains a (^-singular part. 
This fact can be displayed in the explicit way as follows: 

P t (x) = e-^-^5[x - (/ - g)t] + P t {x), (3.18) 

where 

Pt(x) = ±- [°° Mv)e- lvx d V (3.19) 

Z7T J -00 

and 

Uv) = Mv) - e-V-M*- 4 "* (3.20) 

with <pt{v) as \rj\ — > 00. The formulas (I3.17I) - (I3.20I) . which give the probability density 
Pt(x) in a much simpler and transparent form as compared to ( 12.101) and ( 12. lip , constitute 
our main results. To the best of our knowledge, this is the first example of nontrivial system 
with quenched disorder for which the statistical properties of the particle coordinate X t can 
be studied in full detail. 

The above derived findings can be used also for studying the probability density in the 
case of the dichotomous random function g[x) possessing statistical properties different from 
the especially chosen initial condition. In particular, if g(+0) = —g with a probability p± 
and g(+0) = +g with a probability P2 which satisfy the condition pi + p 2 = 1, then the 
corresponding probability density Rt(x) is given by Rt(x) = p\P t (x) + p2Pt{x)\ g ^- g . 



IV. MOMENTS OF THE PROBABILITY DENSITY 

The moments of the probability density, (X^) = §° 00 x k P t (x)dx (k = 1,2,...), are ex- 
pressed through the characteristic function in the commonly known way: 

In principle, with the help of Eqs. (14. ip and (13.171) any moment can be calculated straight- 
forwardly. Specifically, for the first moment these equations yield 

(X t ) = t - ^=^(1 - e- 2 ^). (4.2) 
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In the case of small times, obeying t <C 1/A/, this expression reduces to (X t ) — (f — g)t. 
It has a simple interpretation: At small times the condition g{X t ) = —g holds with almost 
unit probability and so almost with probability one a particle moves with the dimensionless 
velocity / — g. If t ^> g/Xf 2 , then (X t ) = (f 2 — g 2 )t/ f, i.e., the long-time average velocity 
of a particle, (v) = \im t ^ 00 (X t ) /t, equals (f 2 — g 2 )/ f (see Fig. 2, a). In order to elucidate 
this result further, we first note that the total average length of the odd intervals S2m-i (and 
the even intervals S2 m ) on the interval (0, (X t )) tends to (X t )/2 as t — > oo. Accordingly, the 
total time which a particle spends in the odd and even intervals are i_ = (X t )/2(f — g) and 
t+ = (Xt)/2(f + g), respectively. Therefore, using the condition £_ + 1 + = t, one obtains in 
the long-time limit that (v) = (f 2 — g 2 )/ f. We emphasize that this result does not depend 
on the explicit form of the probability density p(s) of the random intervals Sj. 
According to Eqs. (14.11) and (|3.17p . the second moment is given by 

(Xt) - p t — 1 

g 2 (f-g)(f + 3g) + g(f-g) 



2A 2 / 4 2A 2 / 4 

x[2A/(/ 2 + /)t + ^(/ + 3^)] e - 2 ^. (4.3) 

As a consequence, the variance a 2 = (X 2 ) — (X t ) 2 of the particle coordinate X t can be 
represented in the form 



- g2{ / x2f f ] [4A/(/ + g) t - 3f - 5g + 4(/ + g 

-2Xfgt)e- 2 ^-{f-g)e~^ t }. (4.4) 



At small times, when t <C 1/A/, this position variance of grows ast 3 ,of = (A/3)Xg 2 (f —g)t 3 , 
and at large times, when t ^> 1/A/, it grows as t, a 2 = g 2 (f 2 — g 2 )t/Xf 3 . The last 
result evidences that the particles exhibit normal biased diffusion with an effective diffusion 
coefficient _D e fr, reading: 

^,,A^ 2 -/>. (4.5) 
cfi t~oo 2t 2Xf 3 V ; 

Interestingly, D c r varies non-monotonically with the external force /: D c r = (/ — g)/X for 
/ » g, D cS = g 2 /2Xf for / > g, and maxl^ = s/(3 3 / 2 A) for f = VSg (see Fig. 2, b). 
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V. TIME EVOLUTION OF THE PROBABILITY DENSITY 



According to Eqs. (13.201) and (13.171) . at t = the relation ^o(^) = holds and thus 
Po(x) = and Pq(x) = 6(x). This result is a direct consequence of the initial condition 
X = 0. At t > the probability density Pt{x) contains both a 5-singular part and a regular 
part. The 5-singular part arises from the existence of a finite, nonzero probability that the 
random function g(x) does not change at all the sign on the interval [0, (/ — g)t]. The total 
probability of these sample paths, Wo(£) = e~ x ^~ 9 ^, defines the weight of the (^-singular 
distribution part and, because the particles move with the velocity f — g, this distribution 
moves towards the right with the same velocity. 

Because the minimal and maximal velocities of the particles are f — g and / + g, re- 
spectively, the regular part of the probability density, Pt{x), is concentrated on the interval 
[(/ — fi'K (/ + 9)A- If x belongs to this interval, then lim^ Pt{ x ) — M/ — fl0/2<7- For small 
but finite times Pt{x) is an almost linear, decreasing function of x which is transformed to 
a unimodal probability density as t increases. These features of the probability density are 
illustrated in Fig. 3. In order to verify and test the theory for calculating Pt{x), we per- 
formed as well direct numerical simulation of Eq. (12.11) . As depicted in Fig. 4, our analytical 
results are in excellent agreement with the numerical findings. 



A. Asymptotic Approach to a Gaussian Density 



Based on the central limit theorem of probability theory, see, e.g., Ref. 24|, we may 
expect that the probability density Pt(x) tends to the normal (Gaussian) probability density 
as t — > oo. To gain more insight into the long-time behavior of Pt{x), we consider the scaled 
probability density defined by 

VM) = a t P t ({X t ) + a£) (5.1) 
which, in accordance with (13.131) . can be written in the form 

Pt(0 = 7T - f°° e'^^-^du. (5.2) 



2tt 
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(f>t(u/a t ) e - iu{Xt),(7t = e~ u2/2 



(5.3) 



1 /3 

Assuming that \u\ <C cr t and keeping only the first two terms of the asymptotic expansion 
of the function <p t {u/a t ) e - tu ( x t)/°~t as t — > oo, we obtain 

I " 

2^J\f(f 2 -g 2 )t- 

Substituting this formula into (I5.2p yields Vt(£) = Voo(0 + ^t^(O) where 

-Poo(0 = r- / e- u l^<du = -=, (5.4) 



is the probability density of the standard normal distribution and 

.2 /-oo 



Poo(0 (5.5) 



2 V*/(/ 2 - 2 2 )t 

(l^t (OI/^ooCO *C 1) describes the deviation of "Pt(£) from Voo{£)- According to this result, 
Vt{£) slowly approaches (oc t~ 1//2 ) the asymptotic Gaussian form as t — > oo. 

A quantitative measure of non-Gaussian behavior of Pt(x) can be characterized by the 
kurtosis defined as 

k(t)= «*t- W) 4 > _ 3 , (5. 6) 

Notably, fc(t) is equal to zero if X t has the Gaussian distribution. Its dependence on time 
is calculated with the help of Eq. (14.11) and is depicted in Fig. 5. In accordance with the 
central limit theorem, the kurtosis tends to zero as t increases. At very small times the 
formula (I3.18P yields 

P t (x) = [1 - A(/ - g)t] 6[x - (/ - g)t] + ^f^-, (5.7) 

where x G [(/ — g)t, (f+g)t). For this probability density we readily obtain {(X t — (A 4 )) 4 ) = 
(16/5)As 4 (/ - 9)t\ a 2 = (A/3)Xg 2 (f - g)t\ and so k(t) = (9/5) [A(/ - g)t}- 1 as t - 0. The 
divergence of kurtosis as t — > corroborates the fact that at small times the probability 
density indeed strongly differs from the normal one. 

The fact that Pt(x) approaches asymptotically the normal probability density function 
with of oc t as t — > oo leads to an interesting conclusion: Under certain conditions the 
long-time behavior of the considered particles is the same as the long-time behavior of 
Brownian particles. In order to elucidate this, we use the dimensionless equation of motion 

11 



for overdamped Brownian particles: Xt = F+£(t) (xq = 0). Here, Xt is the Brownian particle 
coordinate, F is an external force, and £(t) is Gaussian, thermal white noise with zero mean 



and correlation function £(£)£(£') = 2A5(t — t') (the overbar denotes an average over all 



realizations of and A is the noise intensity). As is well-known, see, e.g., Refs. 



25 



26|, 



Xt has a Gaussian density with xt = Ft and (xt — Xt) 2 = 2At. Therefore, the long-time 
statistical properties of usual Brownian particles and particles in our case are the same if 



x t = (X t ) and (x t — x t ) 2 = of, i.e., F = (f 2 — g 2 )/ f and A = D e g. We further note that in 
the Brownian case the internal noise intensity A is proportional to the absolute temperature; 
thus, the long-time behavior of X t can also be characterized by the dimensionless effective 
temperature T eS = D cS . 



B. Numerical Simulations 



The numerical simulation of Eq. ( 12.11) turned out to be a valuable tool in verifying our 
analytical findings. We take advantage of the fact that this equation contains no explicit time 
dependent terms. In such a case, the numerical simulations can be made especially efficient. 
According to Eq. (12.11) . inside each interval where g(x) = —1 or +1 the particle velocity is 
given by / — g or / + g, respectively. Thus, the algorithm consists of successive generations 
of random interval lengths according to the exponential distribution and a calculation of the 
times needed to travel the generated interval with constant velocity / — g or / + g. The total 
time is summed up until it reaches the required final value. The ensemble average is then 
obtained by repeating this outlined procedure with different realizations of random intervals. 
The probability density function is then obtained as the histogram of final positions of the 
particle. Using this method, we could obtain the probability density Pt(x) with similar 
computational effort as the calculation of the Fourier integral (13. 13j) . The averaging over 
10 7 realizations usually assumed several seconds on modern workstations. The existence of 
the 5-singular part of Pt(x) was recognized in our histogram-procedure by changing the bin 
size of the histogram: The number of counts in the histogram bin containing the 5-function 
does not depend on the bin size. 
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VI. DISCUSSION AND CONCLUSIONS 



We used a path integral approach for determining the probability density function of 
overdamped particles in a piecewise linear random potential which are in addition driven by 
a constant force. Assuming that the intervals of the piecewise linear parts of the potential are 
distributed with an exponential distribution, we succeeded in obtaining the time- dependent 
probability density in the form of an explicit Fourier integral. Within this framework, we 
showed that the probability density contains, apart from a regular part, also a 5-singular 
contribution. The weight of this 5-singular part exponentially decreases with time and the 
total probability density slowly, as t -1 / 2 , converges to a Gaussian density in the long-time 
limit. 

We further calculated the first and the second moments of the probability density function 
and studied its time evolution, both analytically and numerically. Our analytical results are 
in perfect agreement with the numerical ones obtained from simulations of the equation 
of motion. Moreover, we derived a simple representation of the probability density for 
small and large times and, to characterize the non-Gaussianity of the probability density, 
we calculated the kurtosis as a function of time. We showed that under certain conditions 
the long-time behavior of the considered particles is the same as Brownian particles. The 
corresponding effective diffusion coefficient and effective temperature are calculated. 
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APPENDIX A: PROOF OF THE NORMALIZATION CONDITION FOR W n (t) 

Let us introduce the quantities 




(Al) 
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with n > 1. According to the definitions of £l n (t) and Sk, these quantities can be written in 
the form 

S n (t) = p(si)ds 1 p(s 2 )ds 2 ... p(s n )ds n . (A2) 
Jo Jo Jo 

Using in Eq. i \2.8h the integral relation p(s)ds = 1 — J Sn+1 p(s n+ i)ds n+ i resulting from 
the normalization condition for p(s), we can express the probabilities W n (t) through the 
quantities S n (t) as follows: 

W n (t) = S n (t)-S n+1 (t). (A3) 

Taking also into account that S^t) = 0, this representation of W n (t) yields Yl™=i W n (t) = 
S\(t). On the other hand, since Wo(t) = 1 — Si(t), we find that the normalization condition 
W (t) + En=i W n (t) = 1 holds true. 

APPENDIX B: DERIVATION OF EQ. (ETTTh 

According to (13.141) . the characteristic function 4>t{v) can be written in the form 

fafa) = _ e -( A -^)/* [A(r])e {x - iv)9t + B( V )e~ {x - iri)9t ] (Bl) 

with 
and 

B( V ) = — r , A(/ ~ P) r e^Hz. (B3) 
KU 27T J -oo (z - z 3 ){z - z 4 ) K ' 

Applying the residue theorem for calculating the integrals in (1B2I) and (1B3I) . we obtain 

A{rj) = i[Res #(77; z x ) + Res #(77; z 2 )\ 

B(r}) = i[Res $(77; z 3 ) + Res $(77; z 4 )] . (B4) 
Here, Res\l/(r/; z lj2 ) and Res $(77; z 3i4 ) are the residues of the functions 

^; g ) = fe , + y + y A "^ e«*»* 
(z-z 1 )(z-z 2 ) 
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of the complex variable z at the points Zk- Assuming that y > X(f + g) (in this case 
Imzk > 0) and taking into account that Res ^(77; Zip) = lim 2 ^ 2l 2 (z — 2^2)^(77; 2) and 
Res $(77; z 3A ) = \im z ^ Z3A (z - z 3A )®(rj; z), from Eqs. (fB4l) . (|B5|) and f)3.15j) we find 

= - e -( A ~^ 9 * f ~~~~~~~~~r sinh [6(n)t - ia(n)t] 
\o(r/) - ia(r]) 



+ cosh [&(r/)t — ia{rj)t] 

B(rj) = ,, X [ 9 ~ f } e ^ )9t sinh [bMt - iof^tl. (B6) 
0(1]) — iayq) 

Finally, substituting (IB6[) into (IBip . we end up with (13.171) . 
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FIG. 1: Model potential and force field, (a) Schematic representation of the piecewise varying linear 
random potential U(x) and (b) the corresponding dichotomous random force g(x) = —dll(x)/dx 
as functions of the particle coordinate x. 
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FIG. 2: (a) Behavior of the average velocity (v) and (b) of the effective diffusion coefficient D c g 
versus the external force / for the chosen parameter values A = g = 1. 
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FIG. 3: Time evolution of the probability density. Depicted are the plots of the probability density 
P t {x) derived from Eqs. (l3TT7l) - (l33ni ) for the evolution time (a) t = 0.75, (b) t = 2, and (c) t = 5. 
The values of the other parameters are set at / = 1, g = 0.3, and A = 1. The vertical arrows depict 
the singular parts of decreasing weight of Pt(x). 
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FIG. 4: Theory versus numerics. We show the probability density Pt(x) as a function of the 
coordinate x for / = 1, g = 0.3, A = 1, and time t = 4. The solid line (red online) represents 
the theoretical result obtained from Eqs. (I3.17| )-( |3.20p . while the histogram refers to our numerical 
simulation of Eq. (|2.ip . 
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FIG. 5: Deviation from Gaussian behavior. The kurtosis k(t) of the probability density is depicted 
on a linear [k(t)] versus logarithmic [t] scale as a function of the evolution time t for / = 1, g = 0.3, 
and A = 1. 
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